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We study equilibrium sequences of close binary systems on circular orbits and composed of differ- 
ent mass stars with polytropic equation of state in Newtonian gravity. The solving method is a mul- 
tidomain spectral method which we have recently developed. The computations are performed for 
both cases of synchronized and irrotational binary systems with adiabatic indices 7 = 3, 2.5, 2.25, 2 
and 1.8, and for three mass ratios: Mi /A/2 = 0.5,0.2 and 0.1. It is found that the equilibrium 
sequences always terminate at a mass shedding limit (appearance of a cusp on the surface of the 
less massive star). For synchronized binaries, this contrasts with the equal mass case, where the 
sequences terminate instead by a contact configuration. Regarding the turning point of the total 
angular momentum (or total energy) along a sequence, we find that it is difficult to get it for small 
mass ratios. Indeed, we do not find any turning points for M1/M2 < 0.5 in the irrotational case. 
However, in the synchronized case, it becomes easier again to find one for mass ratios much smaller 
than M1/M2 ~ 0.2. 

PACS number(s): 04.25.Dm, 04.40.Dg, 95.30.Lz, 97.80.-d 



I. INTRODUCTION 

Among the binary pulsars observed until now in the Galaxy [QJ^], there exist three double neutron star systems 
which are expected to merge within the age of the Universe (^-Q] . Such close binary neutron stars are important not 
only as a suitable field for testing the predictions of general relativity ^ but also as one of the promising sources of 
grav itational waves for ground based laser interferometers || and as one of candidates of gamma-ray burst sources 

Binary systems of neutron stars evolve due to the emission of gravitational waves, which makes the system lose 
energy and shrinks the orbits. The evolutionary sequence can be separated into three stages. The first one is the 
inspiraling stage in which the orbital separation is much larger than the radius of a neutron star, and the point mass 
post-Newtonian treatment constitutes an excellent approximation fll|| . The second one is the intermediate stage in 
which the orbital separation is only a few times larger than the radius of a neutron star, so that hydrodynamics as 
well as general relativity plays an important role. In this stage, since the shrinking time of the orbital radius due 
to the emission of gravitational waves is still larger than the orbital period, it is possible to approximate the state 
as quasiequilibrium (see fl^-|l5|| for the formulation and []l6|-[20f f or numerical computations). The final stage is the 
merging stage in which the two stars coalesce dynamically |21|-|28(|. 

The present article is the third of a series ^,^9| devoted to the study of the intermediate stage, i.e. the quasiequi- 
librium one, in order to investigate the physical processes around the innermost stable circular orbit (ISCO) as well 
as to prepare initial conditions for computing the merging stage. In Paper I pfj[ (general relativistic computations 
with 7 = 2) and Paper II |p9| (Newtonian computations with various 7), we considered binary systems composed of 
identical neutron stars for calculational simplicity. However, the masses of neutron stars which have been measured 
by the pulsar timing technique including relativistic effects are not exactly the same but span some range jl| . 

In a recent work Jl9[ , Uryu, Shibata, & Eriguchi have shown that the coalescence of irrotational binary neutron 
stars F] of equal mass will produce a disk of negligible mass fraction around black hole. This makes this process 
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It has been shown that the gravitational-radiation driven evolution is too rapid for the viscous force to synchronize the 
spin of each star with the orbit. Rather the viscosity is negligible and the fluid velocity circulation is conserved in these 
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unlikely to be a gamma-ray burst source. As suggested by them, however, the situation may change if one considers 
binary neutron stars of unequal mass. 

Until now, all computations of irrotational binaries of fluid stars have concerned equal mass systems Jl6|- p0| , p2| , p9t 
except for the Newtonian semi-analytic solutions of the ellipsoidal approximation |]33| and the perturbative approach 
pi| , [35| 0. Regarding synchronized binaries, computations of different mass systems have been performed in Newtonian 
gravity |3tj|-|3q| and in the first post-Newtonian regime |39| . In the present paper, we present numerical results of 
synchronized and irrotational binary systems composed of different mass stars with the same polytropic equation 
of state in Newtonian gravity. The solving method is a multidomain spectral method which we have developed in 
both general relativity and Newtonian gravity [ ^0p^j20| , |29| . The results of relativistic calculation will be given in a 
forthcoming paper |41[ |. 

The plan of the article is as follows. We give a short overview of our method in Sec. II. The numerical results are 
presented and discussed in Sec. III. Section IV is then devoted to the summary. 



II. METHOD 



A. Equations to be solved 

The reader is referred to Sec. II. A of Paper II |^9| for the complete set of equations governing perfect fluid Newtonian 
binary stars. Here let us simply mention that, although we consider binary systems composed of different mass stars, 
we assume that the two stars obey to the same polytropic equation of state. This assumption could easily be relaxed 
in our numerical code, in order to study other astrophysical objects, such as neutron star/white dwarf binary systems. 



B. Numerical procedure 



We are using a multidomain spectral method which has been presented in great details in Paper I p0|. For stiff 
equation of states (adiabatic index 7 > 2), we employ the regularization technique described in Paper II p9| in order 
to get rid of the infinite gradient of the density field at the surface of the star and obtain highly accurate results. 

The solution is obtained by means of an iterative procedure, the various steps of which are given in Sec. II. B of 
Paper II. The modifications of this procedure to take into account that the two stars are different concern steps (b) 
and (c), which become 

(b) The separation between the centers of the two stars, d, is held fixed. Here we define the center as the point of 
the maximum enthalpy (or equivalently maximum density). As an initial step, we set the X coordinates of the 
two stellar centers as 

X<1> > ini = ~Mi+M 2 d > (1) 



where Mi and M 2 denote the gravitational masses (or equivalently baryon masses in Newtonian gravity) of the 
two spherically symmetric initial stars. Note here that the location of the rotation axis A rot , which is initially 
set to 0, will change. This is because during the relaxation of the binary system to the equilibrium figure, the 
stars will be deformed and the positions of the centers of masses of each star will change (see Sec. IV.D.l of 
Paper I for details). 



systems. Provided that the initial spins are not in the millisecond regime, this means that close binary configurations are well 
approximated by zero vorticity states |3C|,pl| . 

2 In Taniguchi & Nakamura's papers [ p4pq ]7 although the derived equations include the different masses case, the corresponding 
results are not shown in tables. However, it is easy to obtain the results for the different masses configurations by multiplying 
some functions of the mass ratio. One can find such functions from the source terms of ordinary differential equations which 
are given in MJ33]. 
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(c) By setting the central values of the gradient of enthalpy to be zero for each star, we calculate the location of the 
rotation axis and the orbital angular velocity f2 (see Sec. IV. D. 2 of Paper I for details). After the determination 
of the rotation axis, the stars are translated parallel to the X-axis in order for the rotation axis to coincide with 
the origin of the coordinate system. In short, we set 

A<l> : new = -X<l>,old — -^rot > (3) 
-X<2> : ncw = A"<2>,old _ X ro t , (4) 
-^rot,new — 0. (5) 

Since this procedure is only the transfer of the coordinate location of the binary system, it introduces no change 
in all the other quantities. 

Numerous tests passed by our numerical codes have been presented in Papers I and II. The only additional test we 
provide here for different mass systems is to evaluate the error in the virial theorem (see Eq. (34) of Paper II) and 
to check that it is small. We systematically list the virial error in all the tables given in this article. For the ^150 
configurations listed in Tables | to VI, it is always below 4 x 10~ 5 , except for 2 configurations, for which it is around 
10" 4 . 



III. RESULTS 



By means of the method presented above, we have computed equilibrium sequences of both synchronized and 
irrotational binary systems on circular orbits in Newtonian gravity. In the present paper, we consider the case of binary 
systems composed of different mass stars but built upon the same polytropic equation of state. We use 3 domains 
(one for the fluid interior) for each star and the following number of spectral coefficients: N r x NgX N v — 33 x 25 x 24 
in each domain. 



A. Equilibrium sequences 

The results for evolutionary sequences (constant-mass sequences) with adiabatic indices 7 = 3, 2 and 1.8 are 



presented in Tables | - Vi in the cases of synchronized and irrotational binaries with mass ratios M1/M2 = 0.5,0.2 
and 0.1. 

In these tables, d denotes the separation between the centers of the two stars. Let us recall that we define the center 
of a star as the point of maximum enthalpy (or equivalently maximum density). On the other hand, d,Q denotes the 
separation between the centers of mass of the two stars. Ro, a\, 02, 03, and ai !0pp are values relative to star 1 (less 
massive star) and denote respectively the radius of a spherical static star of the same mass, the radius parallel to 
x-axis toward the companion star, the radius parallel to y-axis, the radius parallel to z-axis, and the radius parallel 
to x-axis in the direction opposite to the companion star. The (2, y, z) axes are the same as in Fig. 1 of Paper I. In 
particular, the stellar centers are located on the x-axis and the z-axis is perpendicular to the orbital plane. The values 
p c and pco indicate the central density of star 1 and that of a spherical static star of the same mass. The prime " ' " 
denotes values relative to the companion star (massive star). The normalized quantities £l, J, and E are defined by 

J:= (GMfi?,,) 1 / 2 ' (7) 

where f2, J, and E denote respectively the orbital angular velocity, the total angular momentum, and the total energy 
of the system, and po is the averaged density of a spherical static star having the same mass as star 1: 



We give also in Tables | - VI the virial error, which provides a measure of the relative accuracy of the numerical 
computation (see Eq. (Ill) of Paper I or Eq. (34) of Paper II for the definition). It has been shown by comparison 
with exact solutions (Roche ellipsoids), that the virial error indicator is very well correlated with the numerical error. 
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Also listed in Tables - VI is the ratio 



X ■■= 



(dH/dr) 



cq.comp 



(dH/8r) polc 



(10) 



where (dH/dr) cq comp [resp. (dH/dr) po \ e ] stands for the radial derivative of the enthalpy at the point on the stellar 
surface located in the orbital plane and looking toward the companion star [resp. at the intersection between the 
surface and the axis perpendicular to the orbital plane and going through the stellar center [z axis)]. This quantity is 
useful because the mass-shedding limit ("Roche limit") corresponds to \ — (cf- Sec. IV. E of Paper I). When \ = 0, 
an angular point (cusp) appears at the equator of the star in the direction to the companion. 

The total angular momentum along an evolutionary sequence is presented in Figs, l] - for the mass ratios 
M1/M2 = 0.5, 0.2 and 0.1. The left panel in each figure is for synchronized binaries, and the right one for irrotational 
binaries. The total angular momenta of irrotational systems have almost the same values as those of point mass 
Keplerian motions because the spin angular momentum of each star is tiny. On the other hand, the angular momenta 
of synchronized systems are larger than those of point mass Keplerian motions because corotating stars have large 
spin angular momentum. 



M,/M 2 = 0.5 

Synchronized case 



Irrotational case 



3.5 



2.5 




FIG. 1. Total angular momentum along an evolutionary sequence for the mass ratio M1/M2 = 0.5. The left panel is for 
synchronized binaries, and the right one for irrotational binaries. Solid, dotted, dashed, long-dashed, and dot-dashed lines 
denote the cases of 7 = 3, 2.5, 2.25, 2, and 1.8, respectively. Note that the scales of the two set of axes are different. 
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M/M 2 = 0.1 

Synchronized case 



11 



Irrotational case 



7.5 




FIG. 3. Same as Fig. [j] but for the mass ratio Mi /Ma = 0.1. 

The relative change in central baryon density along an evolutionary sequence is shown in Figs. ^ - ^| for the mass 
ratios Mx/M^ = 0.5, 0.2 and 0.1. One can see from these figures that, for all mass ratios, the central density decreases 
monotonically as the two stars come closer. This behavior agrees with the analytical results by Taniguchi & Nakamura 




5 



.■£? -0.02 



CD 

rr 



-0.04 



-0.06 



-0.08 



M/M 2 = 0.5 

Synchronized case 




-0.005 



-0.01 



-0.015 



Irrotational case 





A y 

/>■- 
r y 


if 
it 

i 




"1 




// 
ill 

•i$ 

'if 

'/I 
if 
Ms 
13 


Y=3 

y 2.5 

7=2.25 

7=2 

y=1.8 



4 

d„/R„ 



FIG. 4. Relative change in the central baryon density of star 1 along an evolutionary sequence for the mass ratio 
Mi/M 2 = 0.5. Note that the scales of the two set of axes are different. 
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FIG. 5. Same as Fig. ^but for the mass ratio Mi /Ma = 0.2. 
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M/M 2 = 0.1 

Synchronized case Irrotational case 




FIG. 6. Same as Fig. ^ but for the mass ratio Mi /Ma = 0.1. 

Isocontours of the baryon density just before the end point of equilibrium sequences with M1/M2 = 0.5 are presented 
in Figs. - [5] for synchronized binaries and Figs. [l0| - |l2| for irrotational binaries. We can clearly see that there appears 
a cusp on the surface of star 1 (less massive starj/whereas the system is still well detached. In the synchronized case, 
this contrasts with equal mass binaries, for which we have found that the equilibrium sequences terminate by the 
contact between the two stars (Paper II). We will discuss this point in the next section. 



Baryon density (z = 0) 



Baryon density (y=0) 







/ Rr 





R, , 



FIG. 7. Isocontour of the baryon density of synchronized binaries with 7 = 3 and M1/M2 = 0.5 when the separation is 
d/(Ro + R'o) — 1.584. The plots are cross sections of Z — and Y = planes. The thick solid lines denote the stellar surface. 
The small rough on the stellar surface is an artifact of the graphical software. 
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Baryon density (z = 0) 



Baryon density (y=0) 





Baryon density (z = 0) 



Baryon density (y=0) 




-4 



-2 



x / R n 




-4 



-2 



FIG. 10. Isocontour of the baryon density of irrotational binaries with 7 = 3 and M\jMi = 0.5 when the separation is 
d/(Ro + R' ) — 1.599. The plots are cross sections of Z — and Y = planes. The thick solid lines denote the stellar surface. 
The small rough on the stellar surface is an artifact of the graphical software. 
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FIG. 11. Same as Fig. Fol but for 7 = 2 with the separation d/(R + R' ) = 1.677. 
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Baryon density (z = 0) 



Baryon density (y=0) 






x / R Q x / R Q 

FIG. 12. Same as Fig. |ic] but for 7 = 1.8 with the separation d/(R + R' ) = 1.751. 



B. End points of sequences 



A good indicator of the appearance of a cusp at the stellar surface is the quantity x already introduced in Papers I 
and II and defined by Eq. (10|). Comparing x with d/{ai + Oq), which becomes unity for a contact configuration, we 
can determine whether an equilibrium sequence will end by a mass-shedding configuration or by a contact one. 

In Figs, [l^ - [l^, the quantity x is presented as a function of the normalized separation dj (ai + a'i) for 7 = 3,2 and 
1.8. In each figure, we show synchronized and irrotational configurations with the mass ratios M1/M2 — 1,0.5,0.2 
and 0.1. Another view of the quantity x is presented in Fig. 16. In this figure, we fix the mass ratio to 0.5 and change 
7. By extrapolating the curves toward x — in the figures, we find that all the constant mass sequences terminate by 
a mass-shedding point which corresponds to a detached configuration (d/(ai + a[) > 1), except for the synchronized 
case with M1/M2 — 1. For irrotational binaries with M1/M2 = 1, the detached final configuration seems marginal 
from Figs. [ll| - [l5[ However, as discussed in more details Paper II, it seems that it is really the case, at least for small 
values of 7. 

It would not be necessary to extrapolate the curves of Figs. [l| - HI if we could compute up to x = (cusp). 
However, when a cusp appears the stellar surface is no longer differentiable. Now our multidomain method makes the 
boundary of the inner domain fit with the stellar surface [Q. This procedure is essential to get accurate solutions 
in the irrotational case (see Appendix B of Paper I) , and also necessary for the regularization of the density profile 
when 7 > 2 (see Sec. Ill of Paper II). When fitting the inner domain boundary to the stellar surface, the latter is 
assumed to be differentiable, being described by a finite series of differentiable functions of (#, (p). This explains why 
we cannot treat cuspy surfaces. Furthermore, for very close configurations, just prior to the apparition of the cusp, 
the surface is highly distorted so that there appear unphysical oscillations when using a finite series of differentiable 
functions (Gibbs phenomenon). 
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FIG. 13. Equatorial to polar ratio of the radial derivative of enthalpy x as a function of the separation d normalized by the 
total radius ai + a[ in the case of 7 = 3. The left (resp. right) panel is for synchronized (resp. irrotational) binaries. Solid, 
dotted, dashed, and long-dashed lines denote the mass ratios M1/M2 = 1,0.5,0.2, and 0.1, respectively. 
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FIG. 14. Same as Fig. [jj] but for 7 = 2. 
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FIG. 15. Same as Fig. but for 7 = 1.8. 
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FIG. 16. Equatorial to polar ratio of the radial derivative of enthalpy \ as a function of the separation d normalized by 
the total radius a± + a[. The mass ratio is M1/M2 = 0.5. The left (resp. right) panel is for synchronized (resp. irrotational) 
binaries. Solid, dotted, dashed, long-dashed, and dot-dashed lines denote the cases of 7 = 3, 2.5, 2.25, 2, and 1.8, respectively. 
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C. Turning points of the total angular momentum 



In this section, we discuss the turning point of the total angular momentum (and/or total energy^) along an 
evolutionary sequence. This turning point is interesting because it corresponds to the onset of the dynamical instability 
in the irrotational case which is regarded as a realistic rotational state of coalescing binary neutron stars Q], while 
it signals the secular instability in the synchronized case [Em . 



As one can see from Tables Q — [II and Table I of Paper II, there always appears a turning point in the 7 = 3 case 
for synchronized binaries. On the other hand, we do not find any turning point in the irrotational cases for different 
mass binaries. Of course, there exists the possibility that the turning point is located just before the apparition of 
a cusp for slightly different mass binaries. As mentioned above, it is difficult for our method to accurately calculate 
just before the apparition of a cusp because of the Gibbs phenomenon. Therefore, we will discuss the appearance of 
the turning point only in the synchronized case. 

In Fig. [17[ the quantity x & t the turning point of total angular momentum is shown as a function of the mass 
ratio M\j Mi in the synchronized case with 7 = 3,2.5,2.25 and 2. It appears clearly that the 7 = 3 curve has a 
minimum at M1/M2 ~ 0.2 — 0.3. This means that there always exists a turning point of total angular momentum for 
synchronized binaries with 7 = 3. One can also see that the turning point exists marginally in the 7 = 2.5 case at 
around M1/M2 ~ 0.2 — 0.3. On the other hand, it seems that the line of 7 = 2.25 may reach \ — (mass shedding 
limit) around 0.1 < M1/M2 < 0.5. if we extrapolate it. This means that there does not exist any turning point for 
synchronized binaries with 7 < 2.25 and 0.1 < M] I Mo < 0.5. 

Therefore we can draw two conclusions from Fig. [IT]. The first one is that it becomes difficult to find the turning 
point for small 7. The second one is that the mass ratio for which it is more difficult to obtain any turning point is 
about M1/M2 ~ 0.2 — 0.3. This can be explained as follows. For equal mass binaries, the orbital motions of the two 
stars contribute to the orbital angular momentum, which decreases when the separation decreases. In addition, the 
spins of the two stars contribute to the spin angular momentum, which increases when the separation decreases. The 
increase of the spin angular momentum has two parts. The first one is the increase of orbital angular velocity, and 
the second one comes from the increase of the moment of inertia. Then, the total angular momentum has a minimum 
somewhere. Keeping this fact in mind, let us consider the different mass case, for example Mi/M 2 = 0.5. There still 
exists large contribution from the massive star to the orbital angular momentum and spin angular momentum, while 
the moment of inertia of massive star does not increase so much because its deformation is small. Then, since the less 
massive star has to produce much more moment of inertia in order to have the turning point, the deformation at the 
turning point becomes large and \ becomes small (see Table || as well as Table I of Paper II) . This tendency continues 
up to M1/M2 ~ 0.3. However, when the mass ratio becomes smaller than 0.2, the situation changes. Then the center 
of mass of the binary system is located inside the massive star near the turning point. Hence the contribution from 
the massive star to the orbital angular momentum becomes small. In addition, there is a large contribution from the 
massive star to the spin angular momentum because of the synchronous rotation. Accordingly, the less massive star 
does not need to produce large moment of inertia. Therefore, it becomes easier again for the turning point to appear. 

If we apply this explanation to the case of irrotational binaries, we can understand why there does not appear 
any turning point in the irrotational cases of different mass binaries. First of all, let us recall the case of equal mass 
binaries. The orbital motions of the two stars contribute to the orbital angular momentum which decreases when 
the separation decreases. Additionally, the spins of the two stars contribute to the spin angular momentum which 
increases when the separation decreases. These situations are the same as in the synchronized case. However, the 
spin angular momentum is much smaller than in the synchronized case because of irrotation. Next, we consider the 
different mass case. The large difference from the synchronized case is the contribution from the massive star to the 
spin angular momentum. In the irrotational case, it becomes negligible when the mass difference is large, in which 
case the massive star does not deviate from a spherical body so much. Then, in order to obtain the turning point, the 
less massive star must produce much more moment of inertia, i.e., much more deformation. However, before obtaining 
the turning point, the less massive star reaches the mass-shedding point. 



3 The turning points of the total energy and the total angular momentum always coincide (dE — QdJ) Jl3|Jl4| 
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1 1 1 ■ 1 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

Mass ratio M,/M 2 

FIG. 17. Quantity \ a t the turning points of the total angular momentum (and/or the total energy) as a function of 
the mass ratio M1/M2 in the synchronized case. Solid lines with circles, dotted with triangles, dashed with diamonds, and 
long-dashed with squares denote the cases of 7 = 3, 2.5, 2.25 and 2, respectively. Thick lines with filled symbols are for star 1 
(less massive star) and thin lines with open symbols for star 2 (massive star). 

In Fig. |1| we show the orbital separation at the turning point of total angular momentum and at the end point 
of a constant-mass sequence as a function of the mass ratio M\jM% in the synchronized case with 7 = 3. One can 
see from this figure that the orbital separation (normalized by a radius of a spherical static star with the same mass) 
increases when the mass ratio decreases. It also appears in Fig.|l^ that when the mass ratio decreases, the separation 
between the turning point and the end point decreases first and increases again. This is in accordance with the 
behavior discussed above (difficulty of finding the turning point for intermediate mass ratios, M\jMi ~ 0.2 — 0.3). 



Orbital separation at the turning points and the end points 
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6 1 ■ , ■ , ■ , ■ , ■ 
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FIG. 18. Orbital separation at the turning point of total angular momentum and at the end point of a constant-mass 
sequence as a function of the mass ratio M1/M2 in the synchronized case with 7 = 3. Solid line with filled circles and dotted 
one with open triangles denote the turning point and the end point (cusp point), respectively. 
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IV. SUMMARY 



We have studied Newtonian equilibrium sequences of both synchronized and irrotational binary systems on circular 
orbits and composed of different mass stars. 

It is found that the equilibrium sequences for different mass binaries always terminate by the mass-shedding limit of 
the less massive star (cusp point), regardless of the rotation state (synchronized or irrotational). This mass-shedding 
limit occurs for a detached configuration. This contrasts with the case of equal mass synchronized binaries, for which 
sequences terminate by the contact between the two stars (Paper II). This is due to the fact that the tidal force from 
the massive companion star is larger in the different mass case, so that it is easier for the less massive star to be 
disrupted before contact. 

Regarding the turning point of the total angular momentum (or total energy) along an evolutionary sequence, we 
have found that it is difficult to appear for small mass ratio -M1/M2. In fact, for irrotational binaries, we did not find 
any turning point for M1/M2 < 0.5 in the region x > 0.2. Of course, there remains the possibility to find a turning 
point in the region of < \ < 0-2, but it is difficult for our method to calculate in that region because of the Gibbs 
phenomenon. For synchronized binaries, on the other hand, we have found that the turning point is difficult to get 
for mass ratios Mi/M 2 ~ 0.2 — 0.3, but easier again for smaller mass ratios, M1/M2 < 0.2. 
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TABLE III. 
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X 


















a 2 /a 1 


1 ! 1 
O3/O1 


/ / / 
"l,opp/ Oi 


(Pc ~ Pco)/p'cO 


x' 


12.926 


5.000 


4 


974 


0.08242 


11.450 


7 

-35.974 


= 3 

5.562(-6) 


0.9861 
0.9994 


0.9814 
0.9975 


0.9989 
0.9999 


-1.439 
-5.488 


-3) 
-4) 


0.9689 
0.9959 


10.320 


3.993 


3 


950 


0.1155 


10.543 


-36.051 


5.510(-6) 


0.9718 
0.9989 


0.9628 
0.9951 


0.9973 
0.9998 


-2.957 
-1.081 


-3) 
-3) 


0.9376 
0.9918 


9.018 


3.489 


3 


431 


0.1415 


10.107 


-36.099 


5.466(-6) 


0.9566 
0.9983 


0.9435 
0.9926 


0.9950 
0.9997 


-4.644 
-1.626 


-3) 
-3) 


0.9049 
0.9877 


7.815 


3.024 


2 


942 


0.1754 


9.741 


-36.148 


5.402(-6) 


0.9304 
0.9973 


0.9113 
0.9886 


0.9904 
0.9995 


-7.708 
-2.512 


-3) 
-3) 


0.8497 
0.9809 


6.511 


2.520 


2 


388 


0.2310 


9.442 


-36.200 


5.276(-6) 


0.8669 
0.9950 


0.8373 
0.9800 


0.9750 
0.9989 


-1.589 
-4.398 


-2) 
-3) 


0.7207 
0.9663 


f5.776 


2.237 


2 


037 


0.2771 


9.376 


-36.214 


5.158(-6) 


0.7805 
0.9925 


0.7440 
0.9708 


0.9436 
0.9981 


-2.839 
-6.395 


-2) 
-3) 


0.5491 
0.9506 


5.397 


2.093 


1 


801 


0.3078 


9.407 


-36.206 


6.572(-6) 


0.6769 
0.9904 


0.6394 
0.9637 


0.8835 
0.9974 


-4.404 
-7.957 


-2) 
-3) 


0.3379 
0.9382 


5.289 


2.054 


1 


687 


0.3179 


9.434 


-36.201 


2.323(-5) 


0.6087 
0.9896 


0.5733 
0.9611 


0.8218 
0.9971 


-5.190 
-8.510 


-2) 
-3) 


0.1899 
0.9337 


10.014 


5.007 


4 


949 


0.1209 


9.844 


7 

-50.983 


= 2 
4.219(-13) 


0.9751 
0.9998 


0.9672 
0.9989 


0.9973 
1.0000 


-5.466 
-5.090 


-3) 
-4) 


0.9236 
0.9975 


8.021 


4.011 


3 


914 


0.1686 


8.964 


-51.092 


7.547(-13) 


0.9491 
0.9995 


0.9343 
0.9979 


0.9928 
0.9999 


-1.143 
-9.922 


-2) 
-4) 


0.8468 
0.9951 


7.000 


3.500 


3 


364 


0.2069 


8.501 


-51.167 


1.083(-13) 


0.9189 
0.9993 


0.8977 
0.9968 


0.9859 
0.9999 


-1.859 
-1.496 


-2) 
-3) 


0.7616 
0.9927 


5.978 


2.990 


2 


771 


0.2624 


8.046 


-51.259 


1.180(-10) 


0.8525 
0.9988 


0.8220 
0.9949 


0.9648 
0.9998 


-3.492 
-2.412 


-2) 
-3) 


0.5836 
0.9881 


5.588 


2.795 


2 


510 


0.2907 


7.882 


-51.298 


4.957(-9) 


0.8001 
0.9985 


0.7657 
0.9937 


0.9412 
0.9997 


-4.772 
-2.962 


-2) 
-3) 


0.4475 
0.9853 


5.343 


2.674 


2 


308 


0.3112 


7.786 


-51.323 


9.742(-8) 


0.7422 
0.9983 


0.7061 
0.9927 


0.9055 
0.9997 


-6.047 
-3.399 


-2) 
-3) 


0.2982 
0.9831 


5.244 


2.625 


2 


198 


0.3202 


7.749 


-51.333 


4.620(-7) 


0.7021 
0.9982 


0.6662 
0.9923 


0.8735 
0.9997 


-6.753 
-3.601 


-2) 
-3) 


0.2108 
0.9821 


5.210 


2.608 


2 


145 


0.3235 


7.737 


-51.336 


8.403(-7) 


0.6802 
0.9981 


0.6448 
0.9922 


0.8528 
0.9996 


-7.057 
-3.676 


-2) 
-3) 


0.1531 
0.9817 


8.599 


5.000 


4 


898 


0.1519 


9.032 


7 = 
-65.879 


= 1.8 

1.971(-9) 


0.9629 
0.9999 


0.9517 
0.9994 


0.9950 
1.0000 


-1.048 
-3.540 


-2) 
-4) 


0.8728 
0.9985 


6.868 


3.994 


3 


814 


0.2129 


8.173 


-66.013 


1.962(-9) 


0.9204 
0.9997 


0.8998 
0.9989 


0.9853 
1.0000 


-2.278 
-6.959 


-2) 
-4) 


0.7380 
0.9970 


6.030 


3.507 


3 


242 


0.2589 


7.739 


-66.101 


1.903(-9) 


0.8696 
0.9996 


0.8414 
0.9983 


0.9686 
1.0000 


-3.759 
-1.030 


-2) 
-3) 


0.5869 
0.9956 


5.862 


3.409 


3 


117 


0.2702 


7.651 


-66.121 


1.753(-9) 


0.8531 
0.9996 


0.8232 
0.9982 


0.9619 
0.9999 


-4.228 
-1.122 


-2) 
-3) 


0.5397 
0.9952 


5.582 


3.247 


2 


891 


0.2909 


7.506 


-66.156 


1.359(-9) 


0.8159 


0.7831 

D QQ7Q 


0.9440 


-5.245 
-1.301 


-2) 
-3) 


0.4354 

f) QQ44 


5.302 


3.085 


2 


610 


0.3145 


7.362 


-66.194 


1.360(-7) 


0.7511 
0.9994 


0.7162 
0.9975 


0.9018 
0.9999 


-6.743 
-1.521 


-2) 
-3) 


0.2589 
0.9934 


5.246 


3.052 


2 


535 


0.3197 


7.334 


-66.202 


1.308(-7) 


0.7294 
0.9994 


0.6945 
0.9974 


0.8843 
0.9999 


-7.140 
-1.572 


-2) 
-3) 


0.2193 
0.9932 


5.203 


3.028 


2 


463 


0.3236 


7.313 


-66.207 


1.021(-7) 


0.7064 
0.9994 


0.6718 
0.9974 


0.8630 
0.9999 


-7.466 
-1.611 


-2) 
-3) 


0.1599 
0.9930 
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TABLE IV. 



Irrotational case M1/M2 = 0.5 
dG/Ro d/(R + R' ) d/(ai+a[) fi J E Virial error a 2 /ai a 3 /ai ai, op p/ai {p c -p c o)/pcO X 

















a 2 /a 1 


a 3 /oi 


/ / / 

a l,opp/ a l 


(Pc - p'cQ 


)/PcO 


X 


7.515 


3.497 


3.473 


0.09710 


3.166 


7 

-2.623 


= 3 
5.208(-6) 


0.9852 
0.9944 


0.9854 
0.9944 


0.9981 
0.9992 


-4.380 
-6.207 


-5) 
-6) 


0.9746 
0.9902 


6.413 


2.984 


2.950 


0.1232 


2.925 


-2.646 


5.124(-6) 


0.9754 
0.9907 


0.9758 
0.9908 


0.9963 
0.9984 


-1.174 
-1.647 


-4) 
-5) 


0.9576 
0.9837 


5.411 


2.518 


2.467 


0.1591 


2.689 


-2.675 


5.014(-6) 


0.9567 
0.9837 


0.9580 
0.9841 


0.9921 
0.9968 


-3.470 
-4.769 


-4) 
-5) 


0.9250 
0.9712 


4.308 


2.005 


1.910 


0.2245 


2.406 


-2.721 


4.811(-6) 


0.9015 
0.9641 


0.9066 
0.9656 


0.9755 
0.9909 


-1.643 
-2.116 


-3) 
-4) 


0.8257 
0.9349 


4.006 


1.865 


1.745 


0.2508 


2.325 


-2.738 


4.724(-6) 


0.8676 
0.9528 


0.8756 
0.9552 


0.9624 
0.9869 


-2.846 
-3.507{ 


-3) 
-4) 


0.7611 
0.9132 


3.703 


1.725 


1.562 


0.2832 


2.245 


-2.756 


4.595(-6) 


0.8089 
0.9350 


0.8219 
0.9389 


0.9343 
0.9799 


-5.529 
-6.280 


-3) 
-4) 


0.6394 
0.8778 


3.498 


1.632 


1.408 


0.3101 


2.197 


-2.768 


2.786(-6) 


0.7318 
0.9156 


0.7501 
0.9214 


0.8855 
0.9712 


-9.876 
-1.001 


-3) 
-3) 


0.4528 
0.8366 


3.401 


1.590 


1.298 


0.3250 


2.180 


-2.773 


2.041(-5) 


0.6553 
0.9018 


0.6762 
0.9089 


0.8210 
0.9643 


-1.437 
-1.300 


-2) 
-3) 


0.2303 
0.8050 


7.000 


3.500 


3.478 


0.1080 


3.055 


7 

-2.643 


= 2 

1.061(-12) 


0.9852 
0.9964 


0.9853 
0.9964 


0.9978 
0.9995 


-8.155 
-4.977 


-5) 
-6) 


0.9647 
0.9912 


5.979 


2.990 


2.957 


0.1368 


2.824 


-2.667 


1.429(-12) 


0.9754 
0.9940 


0.9758 
0.9941 


0.9957 
0.9990 


-2.165 
-1.301 


-4) 
-5) 


0.9411 
0.9855 


5.007 


2.504 


2.454 


0.1786 


2.586 


-2.700 


5.191(-13) 


0.9556 
0.9894 


0.9567 
0.9895 


0.9905 
0.9979 


-6.673 
-3.885( 


-4) 
-5) 


0.8930 
0.9741 


3.986 


1.993 


1.899 


0.2520 


2.311 


-2.750 


3.427(-10) 


0.8978 
0.9772 


0.9023 
0.9778 


0.9695 
0.9942 


-3.108 
-1.647 


-3) 
-4) 


0.7476 
0.9432 


3.644 


1.823 


1.692 


0.2887 


2.214 


-2.773 


4.678(-8) 


0.8493 
0.9686 


0.8570 
0.9696 


0.9453 
0.9912 


-6.072 
-2.975 


-3) 
-4) 


0.6143 
0.9209 


3.400 


1.701 


1.512 


0.3212 


2.145 


-2.791 


5.162(-6) 


0.7771 
0.9592 


0.7888 
0.9607 


0.8962 
0.9874 


-1.091 
-4.805 


-2) 
-4) 


0.3791 
0.8957 


3.350 


1.677 


1.463 


0.3286 


2.132 


-2.795 


1.852(-5) 


0.7494 
0.9567 


0.7622 
0.9584 


0.8729 
0.9864 


-1.255 
-5.338 


-2) 
-4) 


0.2773 
0.8890 


3.326 


1.665 


1.432 


0.3324 


2.125 


-2.797 


4.053(-5) 


0.7296 
0.9553 


0.7431 
0.9571 


0.8547 
0.9858 


-1.355 
-5.638 


-2) 
-4) 


0.2014 
0.8853 


6.701 


3.516 


3.493 


0.1153 


2.989 


7 = 
-2.677 


= 1.8 
1.871(-9) 


0.9844 
0.9972 


0.9846 
0.9972 


0.9975 
0.9996 


-1.045 
-3.504 


-4) 
-6) 


0.9580 
0.9924 


5.723 


3.003 


2.970 


0.1461 


2.763 


-2.702 


1.857(-9) 


0.9741 
0.9954 


0.9745 
0.9954 


0.9951 
0.9992 


-2.774 
-9.128 


-4) 
-6) 


0.9299 
0.9874 


4.746 


2.491 


2.437 


0.1935 


2.517 


-2.738 


1.811 (-9) 


0.9515 
0.9916 


0.9527 
0.9917 


0.9885 
0.9983 


-9.087 
-2.874 


-4) 
-5) 


0.8678 
0.9769 


3.824 


2.007 


1.905 


0.2680 


2.263 


-2.788 


5.834(-10) 


0.8913 
0.9829 


0.8957 
0.9832 


0.9638 
0.9958 


-3.888 
-1.109 


-3) 
-4) 


0.6959 
0.9523 


3.489 


1.831 


1.682 


0.3080 


2.165 


-2.813 


5.413(-7) 


0.8327 
0.9765 


0.8402 
0.9771 


0.9291 
0.9936 


-7.716 
-2.006 


-3) 
-4) 


0.5055 
0.9337 


3.348 


1.758 


1.559 


0.3280 


2.124 


-2.824 


1.519(-5) 


0.7783 
0.9727 


0.7883 
0.9735 


0.8856 
0.9922 


-1.087 
-2.636 


-2) 
-4) 


0.2930 
0.9227 


3.334 


1.751 


1.542 


0.3301 


2.120 


-2.825 


2.569(-5) 


0.7688 
0.9723 


0.7791 
0.9731 


0.8768 
0.9920 


-1.131 
-2.714 


-2) 
-4) 


0.2534 
0.9214 


3.323 


1.745 


1.527 


0.3318 


2.116 


-2.826 


3.720(-5) 


0.7596 
0.9720 


0.7702 
0.9727 


0.8679 
0.9919 


-1.170 
-2.777 


-2) 
-4) 


0.2149 
0.9203 
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TABLE V. 



Irrotational case M1/M2 = 0.2 
da/Ro d/(R + R' o ) d/(ai + oi) Q J E Virial error a 2 /ai a 3 /ai ai, opP /ai (p c - Pca)/pcO X 

























a 2 /a 1 


03/ai 


a l,opp/ a l 


{Pc ~ PcC 


)/PcO 


X 


9.508 


3 


995 


3 


971 


0.09649 


6 


295 


-10 


7 
885 


= 3 
5.459(-6) 


0.9820 
0.9981 


0.9822 
0.9981 


0.9982 
0.9997 


-6.669 
-7.198 


-5) 
-7) 


0.9696 
0.9967 


8.374 


3 


519 


3 


487 


0.1167 


5 


908 


-10 


920 


5.427(-6) 


0.9730 
0.9972 


0.9735 
0.9972 


0.9968 
0.9996 


-1.467 
-1.560 


-4) 
-6) 


0.9544 
0.9950 


7.153 


3 


006 


2 


960 


0.1479 


5 


462 


-10 


971 


5.380(-6) 


0.9549 
0.9953 


0.9561 
0.9954 


0.9936 
0.9992 


-3.981 
-4.107 


-4) 
-6) 


0.9237 
0.9917 


5.931 


2 


493 


2 


418 


0.1961 


4 


978 


-11 


043 


5.311(-6) 


0.9140 
0.9913 


0.9177 
0.9915 


0.9844 
0.9981 


-1.380 
-1.327 


-3) 
-5) 


0.8539 
0.9844 


5.232 


2 


199 


2 


090 


0.2370 


4 


683 


-11 


098 


5.255(-6) 


0.8616 
0.9866 


0.8694 
0.9871 


0.9687 
0.9968 


-3.415 
-2.981 


-3) 
-5) 


0.7621 
0.9757 


4.705 


1 


979 


1 


813 


0.2786 


4 


456 


-11 


148 


5.179(-6) 


0.7777 
0.9803 


0.7926 
0.9812 


0.9337 
0.9947 


-8.255 
-6.105 


-3) 
-5) 


0.6022 
0.9636 


4.437 


1 


869 


1 


627 


0.3054 


4 


347 


-11 


176 


3.723(-6) 


0.6847 
0.9752 


0.7054 
0.9764 


0.8765 
0.9928 


-1.529 
-9.337 


-2) 
-5) 


0.3933 
0.9535 


4.362 


1 


840 


1 


551 


0.3138 


4 


321 


-11 


183 


3.632(-6) 


0.6343 
0.9733 


0.6565 
0.9747 


0.8346 
0.9920 


-1.928 
-1.066 


-2) 
-4) 


0.2596 
0.9496 


8.021 


4 


Oil 


3 


975 


0.1245 


5 


782 


-13 


7 
312 


= 2 
7.577(-13) 


0.9754 
0.9990 


0.9757 
0.9990 


0.9968 
0.9999 


-2.289 
-3.480 


-4) 
-7) 


0.9421 
0.9977 


7.000 


3 


500 


3 


451 


0.1528 


5 


402 


-13 


357 


1.708(-14) 


0.9617 
0.9985 


0.9625 
0.9985 


0.9941 
0.9998 


-5.369 
-7.936 


-4) 
-7) 


0.9100 
0.9965 


5.979 


2 


990 


2 


916 


0.1936 


4 


994 


-13 


418 


1.624(-12) 


0.9349 
0.9976 


0.9369 
0.9976 


0.9877 
0.9996 


-1.482 
-2.072 


-3) 
-6) 


0.8468 
0.9942 


5.347 


2 


674 


2 


572 


0.2290 


4 


726 


-13 


467 


1.954(-11) 


0.9028 
0.9966 


0.9066 
0.9966 


0.9780 
0.9994 


-3.147 
-4.116 


-3) 
-6) 


0.7703 
0.9917 


5.006 


2 


504 


2 


378 


0.2529 


4 


575 


-13 


498 


4.254(-10) 


0.8746 
0.9958 


0.8802 
0.9958 


0.9677 
0.9992 


-5.018 
-6.190 


-3) 
-6) 


0.7016 
0.9896 


4.616 


2 


309 


2 


136 


0.2859 


4 


399 


-13 


539 


3.578(-8) 


0.8217 
0.9945 


0.8308 
0.9946 


0.9433 
0.9988 


-9.292 
-1.029 


-3) 
-5) 


0.5645 
0.9864 


4.372 


2 


188 


1 


952 


0.3107 


4 


287 


-13 


568 


1.301(-6) 


0.7601 
0.9934 


0.7726 
0.9935 


0.9046 
0.9985 


-1.474 
-1.458 


-2) 
-5) 


0.3826 
0.9835 


4.263 


2 


134 


1 


838 


0.3229 


4 


238 


-13 


581 


1.304(-5) 


0.7077 
0.9928 


0.7220 
0.9929 


0.8610 
0.9983 


-1.878 
-1.719 


-2) 
-5) 


0.2068 
0.9820 


7.182 


4 


002 


3 


951 


0.1470 


5 


471 


-15 


7 = 
497 


= 1.8 
1.945(-9) 


0.9680 
0.9994 


0.9685 
0.9994 


0.9950 
0.9999 


-4.441 
-1.670 


-4) 
-7) 


0.9148 
0.9984 


6.288 


3 


504 


3 


433 


0.1794 


5 


120 


-15 


547 


1.935(-9) 


0.9503 
0.9991 


0.9514 
0.9991 


0.9909 
0.9999 


-1.029 
-3.727 


-3) 
-7) 


0.8679 
0.9975 


5.394 


3 


006 


2 


897 


0.2259 


4 


745 


-15 


612 


1.917(-9) 


0.9150 
0.9985 


0.9178 
0.9985 


0.9805 
0.9998 


-2.805 
-9.451 


-3) 
-7) 


0.7740 
0.9960 


4.932 


2 


748 


2 


603 


0.2586 


4 


539 


-15 


655 


1.331(-9) 


0.8804 
u.yyou 


0.8851 
u.yyoi 


0.9673 
u.yyy 1 


-5.203 

1 A 
-1.0O4 


-3) 

-O ) 


0.6796 

n act a 7 

u.yy4 ( 


4.499 


2 


508 


2 


296 


0.2970 


4 


340 


-15 


703 


1.190(-7) 


0.8200 
0.9974 


0.8280 
0.9974 


0.9363 
0.9995 


-1.026 
-2.877 


-2) 
-6) 


0.5018 
0.9928 


4.313 


2 


405 


2 


127 


0.3167 


4 


253 


-15 


726 


2.921(-6) 


0.7680 
0.9970 


0.7783 
0.9970 


0.8990 
0.9994 


-1.447 
-3.729 


-2) 
-6) 


0.3232 
0.9917 


4.282 


2 


388 


2 


090 


0.3202 


4 


239 


-15 


730 


5.753(-6) 


0.7539 
0.9969 


0.7647 
0.9969 


0.8869 
0.9994 


-1.542 
-3.934 


-2) 
-6) 


0.2708 
0.9915 


4.251 


2 


370 


2 


046 


0.3237 


4 


224 


-15 


734 


1.294(-5) 


0.7352 
0.9969 


0.7466 
0.9969 


0.8694 
0.9994 


-1.651 
-4.105 


-2) 
-6) 


0.1998 
0.9913 
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TABLE VI. 













Irrotationa 


case 


Mi/Ma = 


= 0.1 












d G /Ro 


d/(Ro + R' ) 


d/(ai + ai) 




J 


E 


Virial error 


a 2 /ai 


03/01 


Ol,opp/oi 


(p c — Pco)/PcO 


X 




















a' 2 /ai 


a 3 /oi 


a l,opp/ a l 


(Pc ~ Pco)/p'cO 


x' 


12.926 


5.000 


4 


981 


0.08241 


10.841 


-35 


7 = 
996 


= 3 
5.587(-6) 


0.9860 
0.9994 


0.9861 
0.9994 


0.9989 
0.9999 


-4.162 
-6.530 


-5) 
-8) 


0.9766 
0.9990 


10.320 


3.993 


3 


960 


0.1155 


9.688 


-36 


093 


5.561(-6) 


0.9716 
0.9989 


0.9720 
0.9989 


0.9973 
0.9998 


-1.671 
-2.540 


-4) 
-7) 


0.9525 
0.9980 


9.018 


3.489 


3 


444 


0.1415 


9.057 


-36 


163 


5.542(-6) 


0.9560 
0.9983 


0.9571 
0.9983 


0.9950 
0.9997 


-3.917 
-5.776 


-4) 
-7) 


0.9266 
0.9970 


7.815 


3.024 


2 


960 


0.1754 


8.434 


-36 


248 


5.518(-6) 


0.9289 
0.9973 


0.9315 
0.9973 


0.9903 
0.9995 


-9.961 
-1.393 


-4) 
-6) 


0.8816 
0.9951 


6.512 


2.520 


2 


413 


0.2309 


7.708 


-36 


375 


5.482(-6) 


0.8618 
0.9950 


0.8694 
0.9951 


0.9744 
0.9989 


-3.594 
-4.379 


-3) 
-6) 


0.7690 
0.9909 


5.706 


2.210 


2 


026 


0.2821 


7.236 


-36 


479 


5.431(-6) 


0.7521 
0.9920 


0.7690 
0.9923 


0.9334 
0.9980 


-1.087 
-1.037 


-2) 
-5) 


0.5695 
0.9854 


5.397 


2.093 


1 


818 


0.3075 


7.063 


-36 


523 


4.325(-6) 


0.6482 
0.9902 


0.6710 
0.9906 


0.8696 
0.9974 


-2.031 
-1.524 


-2) 
-5) 


0.3440 
0.9818 


5.343 


2.074 


1 


762 


0.3124 


7.036 


-36 


531 


1.675(-6) 


0.6132 
0.9898 


0.6367 
0.9902 


0.8400 
0.9972 


-2.362 
-1.641 


-2) 
-5) 


0.2549 
0.9810 


10.014 


5.007 


4 


963 


0.1209 


9.542 


-50 


7 = 
999 


= 2 
3.959(-13) 


0.9751 
0.9998 


0.9755 
0.9998 


0.9973 
1.0000 


-2.407 
-2.281 


-4) 
-8) 


0.9421 
0.9994 


8.021 


4.011 


3 


937 


0.1686 


8.542 


-51 


123 


4.836(-13) 


0.9489 
0.9995 


0.9502 
0.9995 


0.9929 
0.9999 


-9.720 
-8.702 


-4) 
-8) 


0.8816 
0.9989 


7.000 


3.500 


3 


395 


0.2069 


7.982 


-51 


214 


3.847(-14) 


0.9187 
0.9993 


0.9215 
0.9993 


0.9862 
0.9999 


-2.370 
-1.986 


-3) 
-7) 


0.8121 
0.9982 


5.978 


2.990 


2 


816 


0.2624 


7.383 


-51 


335 


4.664(-10) 


0.8526 
0.9988 


0.8597 
0.9988 


0.9657 
0.9998 


-7.186 
-5.205 


-3) 
-7) 


0.6571 
0.9971 


5.588 


2.795 


2 


564 


0.2905 


7.144 


-51 


392 


2.349(-8) 


0.8008 
0.9985 


0.8114 
0.9985 


0.9429 
0.9997 


-1.218 
-7.896 


-2) 
-7) 


0.5264 
0.9964 


5.344 


2.674 


2 


372 


0.3109 


6.993 


-51 


431 


5.703(-7) 


0.7448 
0.9983 


0.7584 
0.9983 


0.9091 
0.9997 


-1.797 
-1.044 


-2) 
-6) 


0.3673 
0.9958 


5.245 


2.625 


2 


270 


0.3198 


6.932 


-51 


448 


3.080(-6) 


0.7073 
0.9982 


0.7223 
0.9982 


0.8800 
0.9997 


-2.150 
-1.175 


-2) 
-6) 


0.2491 
0.9955 


5.221 


2.613 


2 


238 


0.3221 


6.917 


-51 


452 


5.224(-6) 


0.6940 
0.9981 


0.7094 
0.9982 


0.8682 
0.9996 


-2.257 
-1.211 


-2) 
-6) 


0.2053 
0.9954 


8.599 


5.000 


4 


924 


0.1519 


8.843 


-65 


7 = 
891 


1.8 

1.975(-9) 


0.9631 
0.9999 


0.9637 
0.9999 


0.9951 
1.0000 


-6.061 
-7.780 


-4) 
-9) 


0.9029 
0.9997 


6.868 


3.994 


3 


859 


0.2129 


7.905 


-66 


037 


1.969(-9) 


0.9214 
0.9997 


0.9238 
0.9997 


0.9857 
1.0000 


-2.565 
-3.015 


-3) 
-8) 


0.7949 
0.9993 


6.030 


3.507 


3 


304 


0.2589 


7.411 


-66 


138 


1.799(-9) 


0.8725 
0.9996 


0.8777 
0.9996 


0.9702 
1.0000 


-6.246 
-6.629 


-3) 
-8) 


0.6670 
0.9990 


5.583 


3.247 


2 


972 


0.2908 


7.135 


-66 


203 


1.143(-8) 


0.8224 
0.9995 


0.8304 
0.9995 


0.9480 
0.9999 


-1.103 
-1.060 


-2) 
-7) 


0.5280 
0.9987 


5.303 


3.085 


2 


715 


0.3143 


D.9o9 


-66 


250 


6.188(-7) 


0.7639 
0.9994 


0.7747 
0.9994 


0.9110 
0.9999 


-1.671 
-1.460 


-2) 
-7) 


0.3423 
0.9984 


5.246 


3.052 


2 


646 


0.3194 


6.923 


-66 


259 


1.827(-6) 


0.7440 
0.9994 


0.7556 
0.9994 


0.8949 
0.9999 


-1.834 
-1.533 


-2) 
-7) 


0.2723 
0.9984 


5.232 


3.044 


2 


627 


0.3207 


6.914 


-66 


262 


2.441(-6) 


0.7379 
0.9994 


0.7497 
0.9994 


0.8896 
0.9999 


-1.879 
-1.583 


-2) 
-7) 


0.2501 
0.9984 


5.224 


3.039 


2 


614 


0.3215 


6.909 


-66 


263 


2.940(-6) 


0.7339 
0.9994 


0.7458 
0.9994 


0.8860 
0.9999 


-1.908 
-1.598 


-2) 
-7) 


0.2356 
0.9984 
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